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Abstract 

Background: Synthetic genetic transistors are vital for signal amplification and switching in genetic circuits. 
However, it is still problematic to efficiently select the adequate promoters, Ribosome Binding Sides (RBSs) and 
inducer concentrations to construct a genetic transistor with the desired linear amplification or switching in the 
Input/Output (I/O) characteristics for practical applications. 

Results: Three kinds of promoter-RBS libraries, i.e., a constitutive promoter-RBS library, a repressor-regulated 
promoter-RBS library and an activator-regulated promoter-RBS library, are constructed for systematic genetic circuit 
design using the identified kinetic strengths of their promoter-RBS components. 

According to the dynamic model of genetic transistors, a design methodology for genetic transistors via a Genetic 
Algorithm (GA)-based searching algorithm is developed to search for a set of promoter-RBS components and 
adequate concentrations of inducers to achieve the prescribed I/O characteristics of a genetic transistor. 
Furthermore, according to design specifications for different types of genetic transistors, a look-up table is built for 
genetic transistor design, from which we could easily select an adequate set of promoter-RBS components and 
adequate concentrations of external inducers for a specific genetic transistor. 

Conclusion: This systematic design method will reduce the time spent using trial-and-error methods in the 
experimental procedure for a genetic transistor with a desired I/O characteristic. We demonstrate the applicability 
of our design methodology to genetic transistors that have desirable linear amplification or switching by 
employing promoter-RBS library searching. 

Keywords: Genetic transistor, Input/Output (I/O) characteristics, Promoter-RBS library. Systematic design 
methodology. Design specifications 



Background using interchangeable natural components which are 

Synthetic biology aims to perform various specific func- nonexistent in natural biology [1]. A lot of the recent 

tions in organisms by inserting a designed gene network, literature focuses on performing electronic circuit 

In the past, synthetic biology could be classified as hav- behaviors in organisms using genetic devices such as 

ing two broad purposes. The first was to create artificial toggle switches [2-5], oscillators [6-13], pulse generators 

life from natural biology using the synthetic methods. [14,15], logic gates [16-19], and filters [20-22]. Synthetic 

The other was to assemble some functional components biologists also design various types of genetic circuits 

with different functionalities by employing genetic de- 
vices to solve useful tasks, such as biosensor decisions or 



* Correspondence: bschen@ee.nthu.edu.tw 

^Equal contributors edge detection [23,24]. 

Yab of Control and Systems Biology, Department of Electrical Engineering, In many industries, SUch as the electronics and manufac- 

National Tsing Hua University, Hsinchu 30013 Taiwan _ _ industries, the characterization and Standardization 

Lab of Genetic Circuit Design, Department of Electrical Engineering, ^ 

National Tsing Hua University, Hsinchu 3001 3, Taiwan of components and the institution of specifications are 

Full list of author information is available at the end of the article 



O 



© 2013 Lee et al.; licensee BioMed Central Ltd. This is an open access article distributed under the terms of the Creative 
BIoIVIGCI CGntrsI commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



Lee et al. BMC Systems Biology 201 3, 7:1 09 
httpy/www.biomedcentral.com/l 752-0509/7/1 09 



Page 2 of 1 1 



already key elements in the production line. In the synthetic 
biology, the expression of a specific protein needs a pro- 
moter, a ribosome binding side (RBS), a protein coding se- 
quence and a terminator, which are DNA fragments. The 
Registry of Biological Standard Parts (http://www.partsreg- 
istry.org), formed by MIT, shows many standard BioBricks, 
and these standard biological components provide synthetic 
biologists with a quick and standardized way of construct- 
ing gene circuits. Also, BioFAB provides some sorts of bio- 
bricks (see http://biofab.org/data), enabling the rapid design 
and prototyping of genetic constructs. However, in the past, 
the strength of a promoter and an RBS, which are the main 
components of transcription and translation, were defined 
according to their relative strength with other promoters 
and RBSs. Now, however, the promoter-RBS strength can 
be quantified by measuring the fluorescence of proteins 
whose coding gene is constructed at the downstream of the 
promoter-RBS component. 

Based on the kinetic strengths of promoter-RBS compo- 
nents, promoter-RBS libraries are constructed for gene 
circuit design. In our gene circuit design, promoter-RBS 
libraries are built based on kinetic parameters of the dy- 
namic gene regulation, which are identified by the nonlin- 
ear least squares method based on experimental data. We 
formulate the design specifications of desired gene circuits 
in advance and choose an adequate set of promoter-RBS 
components from the promoter-RBS libraries based on the 
characterized, standardized and quantified components. 
For our promoter-RBS libraries, we select three kinds of 
promoters, i.e., constitutive promoters, repressor-regulated 
promoters and activator-regulated promoters, to combine 
with RBSs as the promoter-RBS components. Each of these 
is constructed using the green fluorescent protein in 
Escherichia coli. and characterized to allow for the con- 
struction of the following three types of promoter-RBS 
libraries: constitutive promoter-RBS libraries, repressor- 
regulated promoter-RBS libraries and activator-regulated 
promoter-RBS libraries. 

To date, several researchers have demonstrated that 
synthetic gene circuits have the functionality of amplifi- 
cation or switching [25-29]. These gene circuits can 
amplify the input signal or switch the output signal as 
it exceeds a specific threshold level. Often shown in 
these genetic amplifiers is the use of a two-stage cascade 
of promoters to achieve the function of amplification. 
However, these circuits only amplify a low level of input 
signal or low concentration of inducer at the first stage 
whfle the second stage consists of the promoter-RBS 
activity being fixed. On the other hand, switching 
circuits switch the output signal using an external in- 
ducer. When the inducer is externally increased, the cir- 
cuit is on, and vice versa. Nevertheless, at present, the 
on-state, or high level, of switching has not been clearly 
defined. 



In this paper, we demonstrate that a simple repressive 
gene circuit can work like an electrical transistor as an 
amplifier or a switch. The amplification gains or switch 
levels of the genetic transistor are regulated by the 
concentrations of inducer and different combinations of 
promoters with RBSs. For the convenience of measure- 
ment and application, reporter genes are constructed as 
the measurable input and output. We show that the I/O 
characteristics of the repressive gene circuit regulated by 
inducer concentration can be effectively predicted by 
adequate selection of promoter-RBS components from 
our libraries. Thus, based on these promoter-RBS librar- 
ies, a look-up table is built to quickly select adequate 
promoter-RBS components for the design of genetic 
transistors with different design specifications. 

In the following sections, we first construct the 
promoter-RBS libraries based on the promoter-RBS 
strength through the dynamic regulatory model of 
promoter-RBS components. Then, we describe the I/O 
characteristics of a genetic transistor with different 
kinetic strengths of promoter-RBS components in the 
promoter-RBS libraries and different concentrations of 
inducers. Finally, a look-up table (or genetic transistor 
library) is constructed for genetic transistor design re- 
quiring prescribed I/O characteristics, which is used by 
searching the most appropriate sets of promoter-RBS 
components and concentrations of inducers via the 
genetic algorithm (GA). 

Methods 

Construction of the promoter-RBS libraries for genetic 
transistors 

In this section, we introduce the characterization and 
standardization of promoter-RBS libraries and employ 
a dynamic mathematical model to construct the 
promoter-RBS libraries according to the identified kin- 
etic strengths of promoter-RBS components, populated 
via experimental data. 

Promoter-RBS libraries based on the identified kinetic 
strengths of promoter-RBS components 

In a systematic design procedure, the characterization 
and standardization of components are important pre- 
paratory tasks before practical design process. These can 
save designers a significant amount of time and avoid 
unnecessary trial-and-error attempts. In the field of syn- 
thetic biology, a particular technique was developed to 
create standard interchangeable biological components 
called BioBricks [30,31]. These allow the synthetic biolo- 
gists to focus on the design of more complex genetic 
circuits rather than the basic construction of the gene 
components. 

BioBricks are DNA fragments with specific functions, 
and include promoters, ribosome binding sites (RBS), 
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repressors, activators, reporters and terminators. In the 
database, there are only a few BioBrick components that 
are well-characterized. The well-characterized BioBrick 
components are conducive to the systematic design of 
synthetic genetic circuits. In order to facilitate the 
design of synthetic gene circuits, wider libraries of well- 
characterized BioBricks need to be constructed. 

In our promoter-RBS libraries, the library indexes are 
the kinetic strengths of promoter and RBS, which are 
considered together as a promoter-RBS component be- 
cause the gene expression is regulated by a promoter-RBS 
component. The kinetic strength of a promoter-RBS com- 
ponent can be systematically identified by a stochastic 
model which simulates the dynamic behavior of promoter- 
RBS components under some external molecular or envir- 
onmental noises. In order to identify the kinetic strength of 
a promoter-RBS component, the green fluorescence gene is 
embedded into the downstream of the promoter-RBS 
component. By measuring the fluorescence dynamic time 
profile and using the nonlinear least squares method [32], 
we identify the kinetic strengths of the promoter-RBS com- 
ponents to be used as the indexes of promoter-RBS 
libraries. 

The construction procedure of the promoter-RBS 
libraries can be generally divided into four steps [33]: (i) 
choose the required promoter-RBS components, (ii) se- 
lect the suitable reporter protein and growth conditions, 
(iii) measure the time-profile data of the dynamic behav- 
ior, and (iv) construct the dynamic regulatory model for 
identifying kinetic strengths of promoter-RBS compo- 
nents to be used as library indexes according to the 
nonlinear least squares method. In the first step, some 
promoters can be regulated by specific transcription fac- 
tors, and different combinations of promoters with RBSs 
give different kinetic strengths of promoter-RBS compo- 
nents, which increase the diversity of the libraries. In 
order to rapidly obtain a variety of kinetic strengths of 
promoter-RBS components, a mutation technique was 
used to create different kinetic strengths of promoters 
and RBSs to increase the varieties of promoter-RBS 
components through the mutation of a specific region 
on promoters or RBSs [18,33]. In the second and third 
step, since different reporter proteins, such as the green 
fluorescent protein or red fluorescent protein have dif- 
ferent degradation rates, the measurement times may 
differ. Further, the cell growth conditions have an effect 
on the results of the measurement. Biological compo- 
nent can be characterized at different cellular growth 
phases, under different culture conditions, or at different 
resolutions. In our experiment, the GFP is selected as 
the reporter protein and the time profiles of fluorescence 
are measured by the microplate reader. In the final step, 
a mathematical dynamic model is built to describe the 
time profile of protein expression. Using the protein 



expression time profile measurements, the nonlinear 
least squares method is employed to identify the kinetic 
strengths of promoter-RBS components to be used as 
the library indexes with the mathematical model. For the 
systematic design of genetic transistors, we construct 
three kinds of promoter-RBS libraries, i.e., constitutive, 
repressor-regulated and activator-regulated promoter- 
RBS libraries. The promoter-RBS components in promoter- 
RBS libraries and all BioBrick components used in this 
study are listed in Additional file 1, respectively. The de- 
tailed construction procedures of constitutive, repressor- 
regulated and activator-regulated promoter-RBS libraries 
are described in Additional file 1. 

Construction and design of the genetic transistor 

After the introduction of regulatory functions of 
promoter-RBS components and the construction of the 
promoter-RBS libraries, we design a synthetic genetic 
circuit, similar to a transistor, with prescribed I/O 
characteristics of amplification or switching through 
the external inducer. Before the construction of the 
synthetic genetic transistor, we introduce the simply 
operation of an electronic transistor in Additional file 1 
to which the genetic transistor will be designed accord- 
ingly following. 

Construction of the genetic transistor 

A genetic transistor is shown in Figure 1(a). The transistor 
is constructed to obtain the output protein concentration 
Xprotein the transistor for amplification or switching 
behavior. The genetic transistor consists of the repressor- 
regulated promoter-RBS component C3 and a repressor 
coding gene. The input repressor Xy.epressor2 to the 
genetic transistor is controlled by the repressor-regulated 
promoter-RBS component C2, which is regulated by the 
corresponding repressor. The input repressor Xrepressor2 will 
form the complex and restrict the production of output 
protein Xprotein by binding the corresponding repressor- 
regulated promoter-RBS component C3 to decrease its kin- 
etic strength. However, when the inducer is added, this 
inducer will bind input repressor Xyepressori and prevent it 
from binding to the repressor-regulated promoter- 
RBS component C3. Then, both the kinetic strength of 
repressor-regulated promoter-RBS component C3 and the 
production of output protein Xprotein will increase. The 
dynamic model of a genetic transistor is described as 
follows: 

^repressor! (^2 5 ^) — P repressor {Pm,C2 1 ^m,C2 1^1^^ 

-(^ -\- YrgprQssor2^Xrepressor2{c2, t) 

(1) 

^protein (^3 5 ^) — P repressor (-^^,C3 1 -^m^cs •> ^repressor2 ■> ^2) 
~ (/^ + y protein) ^protein{C'i, t) 
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(a) 



Input device 



Genetic transistor 



Ouput device 




protein 

'prote,n= OUtpUt Signal 



(b) 



Input signal generation device 



Genetic transistor 



Output signal measure device 



measured output ) 




Input signal measure device 

Figure 1 The representation of synthietic genetic transistor circuit, (a) A genetic transistor, (b) A genetic transistor witli measurement 
circuit. Tine input signal of tine genetic transistor is measured by RFP reporter and tine output signal of the genetic transistor is measured by 
GFP reporter. 



where Xy.epressor2 ^nd Xprotein denote the concentrations of in- 
put repressor2 and the output protein of the genetic tran- 
sistor, respectively, and jprotein denotes the degradation rate 
of the protein. 

However, the protein concentration is difficult to dir- 
ectly measure and quantify. To determine characteris- 
tics of the synthetic genetic transistor, a genetic 
transistor with measurement circuit is constructed as 
shown in Figure 1(b). In Figure 1(b), we construct an 
additional repressor-regulated promoter-RBS compo- 
nent c-i so that the input reporter protein Xx can be 
measured by input fluorescence and the output re- 
porter protein x^ can be measured by output fluores- 
cence gi. Note that RFP is used to measure input while 
GFP is used to measure output. Additionally, for the 
convenience of the input regulation, we construct an 
input signal generation device with the concentration 
of inducer /i to control the input of the genetic tran- 
sistor circuit. Then, the dynamic model of a synthetic 
genetic transistor circuit with I/O measure devices 
under environmental disturbances is described by the 
following set of equations: 



^repressor\ (f\ 5 ^) — PcomA^Ci ) + Y repressor!) ^repressor\ (^1 , t) 

^repressor2 (^2 ? ^) P repressor {,-^-^,^2 1 -^^,02 1 ^repressorl 1 ^1 ) 

-if^ + rrepressor2)^repressor2{c2, t) + V2{t) 
(^2 5 i) = P repressor ij^M,C2 1 ^m,C2 1 ^repressorl 1 ^l) 
< -(mi + ^ + r,m,;ci>lfe, t) + V3(0 

g^{c2,t)=mi -Xi (C2 , 0 - (^ + Ym^x, kl (^2 , 0 + ^4 (0 
^2 (^3 5 t) = P repressor {-^^,03 1 Pm,C3 1 ^repressor2 5 ^2) 

-(^2 + + r,m,;c2)^2(c3, t) + Vs(t) 
^2(^3, t) = m2'X2{c3, t)-{fi^ Ym,x^2,2^^'i^ 0 + ^6(0, 
C\^Lihconsti C2 C-^^Librepressor 

(2) 

where mi and m2 denote the maturation rates of 
reporterl Xi and reporter2 X2, respectively, and Vi(^), 
/ = 1, 2, • • 6 denote the noises. 
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To explore the I/O characteristics of a synthetic gen- 
etic transistor with the function of amplification or 
switching, the steady state model of (2) is given by 



Yrepressorl i^l) — P const ci) / {/^ + Yrepressorl) ~^ ^si 
^repressor! {^2 -j^l) — P repressor (-^^,^2 -> ^m,C2 -> Yrepressorl 5 ^1 ) / 

{H" + Yrepressorl) ~^ 
Yl{c2i Ii) P repressor {pM,C2 •> ■^ni,C2 •> Yrepressorl 7 ^1 ) / 

(mi +i^ + r,m,;ci) +^53 
^i(c2,/i) = mi-xi{c2)l{li + r^^^J + Vs, 

Y2 {^3 1^11^2) = P repressor (j^M,C3 -> ^n^jC^ -> Yrepressor2 1^2) / 

Ym,X2} ~^ 

Ci^Lib const 1 ^2 ^rid C^^Lib repressor 

(3) 



translation rates as stochastic processes in the following 
model: 

Pc.^Pc, + APcMt).PM,C2^PM,C2 + APm,C2W2(0, 
P m,C2 ^P m,C2 

+ APm,C2n2{t),PM,c,^PM,c, + APM,c,n3{t), 
P m,C3 ^P m,C3 ~l~ ^P m,C3 ^3(^)5 
Yrepressorl ^Yrepressorl ~^ ^Yrepressorl ^lip) 1 
Yrepressor2 ^ Yrepressor2 ~^ ^ Yrepressor2 ^2(^)5 

Yim,xi ^Yim,xi ~^ ^Yim,xi^'^i^) ^ Ym,xi ^Ym,xi ~^ ^Yim,xi^'^i^) ^ 
Yim,X2 ^Yim,X2 ~^ ^Yim,X2^^^^) "> Ym,X2 ^Ym,X2 ~^ ^Ym,X2^^^^) "> 

mi^mi + Amin2{t), m2^m2 + Am2n3{t), 
pL^fi + Afini{t) 

(5) 

where AP^ , APm,C2 > ^^m,C2 > APm,C3 > ^Pm.cs > ^YrepressorV 
^Yrepressor2> ^Yim,x,> ^Ym,x,> ^Yim,X2> ^Ym,X2> ^^1' and 

A/u denote the standard deviations of stochastic parame- 
ters to be tolerated and could be specified before design 
and Hiit), i = 1, 2, 3 denote Gaussian noises with zero 



AP. 



where v,^ , / = 1,2, ",6 denote the noises at the steady ^f^^ ^^^^ variance. Therefore, AP,^, APm,c2. ^Pn 
state. 

From (3), if - ^2, Ym,xrYm,x2 nm,^i"nm,:.2' th^^ 
the I/O characteristic can be regarded as input/output = 

Yrepressor2lYprotein^YilX2^gllg2> i.e., We COuld USe the xj 

X2 or gilg2 ratio to replace the I/O characteristic of the 
synthetic genetic transistor. Further, the I/O characteris- 
tic can be controlled and regulated by the selection of 
promoter-RBS components C3 and inducer concentration 
/2. Therefore, we need to define the I/O characteristic of 
synthetic genetic transistor circuits to design a genetic 
transistor with the desired I/O characteristic. This is 
done as follows: 



3^55(^3, ^2, ^i(c2, /i)) =^2fe, A, /2),^ie[^ie, gm] (a-u.; 



(4) 



where jssics, gi) denotes the I/O response of the synthetic 
genetic transistor circuit between input signal ^1 and output 
signal g2, and and gi^ denote the lower bound and upper 
bound of ^1. a.u. stands for arbitrary unit. 

In Figure 1(b), promoter-RBS components Ci and C2 
can be selected to control input signals Yrepressorii^i) 
and Yrepressor2{c2> h) 1^ (3). In general, genetic compo- 
nents are inherently uncertain in the biological system 
as a result of gene expression noises in transcription or 
translation processes, thermal fluctuations, DNA mu- 
tations, evolutions, context-dependence between pro- 
moters, 5'UTRs, and coding sequences, as well as 
parameter estimation errors [34-37]. Hence, we model 
the uncertain kinetic strengths of promoter-RBS com- 
ponents, degradation rate of proteins and transcription/ 



M,C3 ) AP;^ C3 f ^Yrepressorl) ^Yrepressor2) ^Yim,xi ' ^Ym,Xi ' 

Ay-^^^, Ay^^^, Ami, Am2 and A// denote the determinis- 
tic parts of parameter variations and ni{t), / = 1, 2, 3 de- 
note different random fluctuation sources. For robust 
design of the genetic transistor circuit, these parameter 
fluctuations in (5) will henceforth be considered in the 
design procedure so that the synthetic genetic transistor 
can tolerate these kinds of parameter fluctuations in vivo. 
With fixed concentration of inducer I2, we expect that 
the input signal ^1 /output signal g2 (I/O) characteristics 
of the synthetic genetic transistor in (4) would be similar 
to the voltage I/O characteristics of the electronic tran- 
sistor shown in Additional file 1. When the inducer 
concentration Ii increases, the kinetic strength of 
promoter-RBS component C2 increases along with the 
fluorescence of the input signal gi, which means that the re- 
pressor concentration Xrepressofi increases. Due to the fixed 
concentration of inducer /2, the redundant repressors 
Yrepressor2y which are not bound by the inducer I2, will re- 
press the promoter-RBS component C3, and the fluores- 
cence of output signal g2 will decrease. Therefore, the I/O 
characteristic of the synthetic genetic transistor is similar to 
Additional flle 1. Additionally, from Additional file 1, we 
see that if input signal is in the operation range of linear 
ampliflcation, the input signal would be inversely amplifled. 

Now, consider the alternative viewpoint, /.a, the volt- 
age I/O characteristics of an electronic transistor. When 
7^2/^1 increases, the reverse amplification gain will be- 
come large and the operation region of linear amplifica- 
tion wiU narrow as shown in (B1)-(B3) and Additional 
file 1. In the synthetic genetic transistor, we expect that 
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when the concentration of inducer I2 changes as per the 
7^2/^1 ratio in (B2)-(B3), the I/O characteristics would be 
similar to the voltage I/O characteristics of electronic 
transistor in Additional file 1. Due to different concen- 
trations of inducer I2, the effect of the inducer on the in- 
put repressor can vary. When the inducer concentration 
I2 decreases, the I/O characteristics would sharpen, so 
the reverse amplification gain becomes large in the oper- 
ation region of linear amplification. 

Finally, when R2/R1 is large enough in (B2)-(B3), the 
operation region of linear amplification will become too 
narrow and result in a sharp change in this region. Cor- 
respondingly with a synthetic genetic transistor, when 
the inducer concentration I2 is low enough, the input 
signal ^1 will produce a small variation, and the output 
signal g2 will have an acute change like a switch. There- 
fore, according to the analysis above, we could obtain 
varying reverse amplification gains and switch levels by 
changing the concentration of inducer /2. 

Systematic design of a genetic transistor based on design 
specification 

According to the above analysis in Figure 1(b), we can 
obtain different reverse amplification gains or switch be- 
haviors via regulation of different concentrations of in- 
ducer /2. Additionally, due to the output signal g2 being 
under the controlled by promoter-RBS component C3, 
we could change the output range by selecting different 
repressor-regulated promoter-RBS components C3 from 
the repressor-regulated promoter-RBS libraries. In this 
way, we can control the I/O characteristics of a synthetic 
genetic transistor to obtain different reverse amplification 
gains or switch levels by choosing different concentrations 
of inducer I2 and selecting different repressor-regulated 
promoter-RBS components C3 from the repressor-regulated 
promoter-RBS libraries. 

In Figure 1(b), the input signal generation device con- 
sists of a constitutive promoter-RBS component Ci, and a 
repressor-regulated promoter-RBS component C2 and an 
inducer /i. The constitutive promoter-RBS component C2 
is selected to produce the input repressor continually. 
Further, for convenience of design, the repressor-regulated 
promoter-RBS component C2 is selected from the corre- 
sponding promoter-RBS library to have sufficient kinetic 
strength to obtain an adequate maximum regulation range 
of input signal regulated by inducer /i. However, the oper- 
ation region of linear amplification is still limited in the 
amplifier design of genetic transistor, and the input signal 
range might not be fully contained in the operation region 
of linear amplification. Therefore, the input signal range 
should be considered in relation to the design purpose. In 
the procedure of amplifier design, the input operation 
range ^1 g [g^ i, g^ J can be set by transforming the inducer 



concentration Ii into the input fluorescence according to 
(2) or (3) as follows 
Input operation range: 

h^[luJi,u]^gi^[gi^i.gi^u] (a-u.) (6) 

where Ii and gi denote the inducer concentration and 
input fluorescence, respectively, Iij and Ii^j^ denote the 
lower and upper bound of inducer concentrations, re- 
spectively, and gij and gi^j^ denote the lower and upper 
bound of input fluorescences respectively. 

Note that, in the future, when the promoter-RBS li- 
braries are large enough, the promoter-RBS components 
Ci and C2 can be designed and selected to match the in- 
put operation range. However, due to the limited size of 
our promoter-RBS libraries and for the convenience of 
design, we will select the repressor-regulated promoter- 
RBS component C2 from the corresponding promoter- 
RBS library. 

From the above analysis, the design purpose of an 
amplifier will lead to the selection of a suitable 
repressor-regulated promoter-RBS component C3 from 
the repressor-regulated promoter-RBS libraries and con- 
centration of inducer I2, Le., {cs.h}) so that the I/O 
characteristics of the synthetic genetic transistor in (4) 
in a specific input range ^1 e \gij-,gi,u\ can match the de- 
sired I/O response similar to (B2), i.e., 

(7) 

where ^1 / and ^2,m denote the lower bound of input 
fluorescence ^1 and upper bound of output fluorescence 
g2, respectively, and gain denotes the amplification gain 
of the genetic transistor. 

On the other hand, the switching behavior will occur 
when the input signal has a small variation (see Additional 
file 1), i.e., a high level signal can be switched into a low 
level signal and vice versa. In the switching behavior of 
synthetic genetic transistor, each promoter-RBS compo- 
nent has its own basal level. Thus, when the input signal 
increases, the output signal will rapidly decrease to the 
basal level. Therefore, the desired I/O response of a 
switch is described as follows: 

yAgi)=Ls^-^^^f^.gAgw giu] (a.u.) (8) 
^^{gilgt) 

where and denote the high level and low level of 
switching, respectively, and gt denotes the transition point 
of input fluorescence. Moreover, the input signal range of 
I/O characteristics of the switch can be set by (6). 

Finally, for matching the desired I/O response of an 
amplifier or switch, the genetic algorithm (GA) is 
employed to select an adequate repressor-regulated 
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promoter-RBS component C3 in the repressor-regulated 
promoter-RBS libraries and the concentration of inducer 
I2 to minimize the following cost function [38], respect- 
ively, i.e., 



min J{c3,h) 

C3&Lib repressor, [hj,h,u\ 

= min E 

C3&Lib repressor [h,l ,h,u] - 



(9) 



To summarize the above design procedure of a bio- 
logical amplifier and switch, a genetic transistor design 
procedure of by the promoter-RBS library searching 
method using GA is proposed as follows [38]: 

1. Construct the genetic transistor circuit such as in 
Figure 1(a). 

2. Build the dynamic and steady state mathematical 
model in (2) and (3), respectively. 

3. Provide the design specification of amplifier with the 
desired I/O response as in (6) and (7) or switch with 
the desired I/O response as in (6) and (8). 

4. Provide the standard deviations of parameter 
fluctuations and environmental disturbances to be 
tolerated in vivo in (5). 

5. Minimize the cost function /(C3, 12) in (9) by 
selecting an optimal set {C3, 12} via GA. 

Based on the design procedure of a genetic transistor 
using the promoter-RBS library searching method with 
GA, the promoter-RBS component C3 is selected from 
the corresponding repressor-regulated promoter-RBS li- 
brary and the inducer concentration I2 is selected within 
V2,hh,u\> while the cost function is calculated in each 
iteration of the selection process. Then, GA would 
select the most adequate promoter-RBS component C3 
from the corresponding repressor-regulated promoter- 
RBS library and inducer concentration I2 e [I2J, l2,u\ to 
minimize the cost function. 

Results 

In silico synthetic genetic transistor design examples 
based on promoter-RBS libraries 

We have presented the construction and design proced- 
ure of a synthetic genetic transistor. In this section, the 
synthetic genetic transistor is designed and simulated to 
verify the I/O characteristics of amplification and switch- 
ing. Subsequently, based on our promoter-RBS libraries, 
the amplification gain in a specific input operation range 
and switching level are designed by employing GA to 
select the most adequate promoter-RBS components 
and inducer concentrations. Finally, to support future 
application of this method, a look-up table for genetic 



transistors is built for different genetic transistor design 
specifications. 

Amplifier design example of synthetic genetic transistor 

Consider the amplifier design of the synthetic genetic 
transistor. Firstly, to obtain the I/O characteristics of 
amplifier, promoter-RBS components {ci, C2} = {/e, ^3} are 
selected to obtain the maximum input operation range. 
The dynamic model and the steady state model have 
been described in (2) and (3). The input operation range 
and desired I/O response of genetic transistor are speci- 
fied as follows: 
Input operation range: 



/lG[l.5*10-^4*10-^] (mM)^^iG[298,431] (a.u.) 



and 



J^fei) = -2-^1 + 1286 



(10) 



(11) 



where -2 is the desired amplification gain as shown in 
Figure 2. Note that the standard deviations of parameter 
fluctuations that are supposed to be tolerated in vivo are 
given by 



^YLacI = ^-^^YLacD TetR = ^-^^YTetR 
^Yim,x, = ^■^^Yim,x,^^Ym,x, = ^■^^Ym,x, 
^Yim,x, = ^■^^Yim,x,^ ^Ym,X2 = ^•^^Ym,x, 

Ami = 0.05mi, Am2 = 0.05^2, A\i = 0.05A|i 



2,3 



(12) 



and the environmental disturbances Vi{t) are independ- 
ent Gaussian noises with zero mean and unit variance. 
Finally, GA is employed to search a set {cs^IaTcl from 
corresponding libraries to minimize the following cost 
function: 



min JA{c3,IaTc) 

cs^Libret, h^l. 5*10-^, 4*10-^] 



min E 

C3^LibTet Ml -5*10-^, 4*10-^] , 



{yss(c3jaTc,gl)-yAgl)f^gl 
(13) 



Then, the most adequate promoter-RBS component 
from the corresponding library and aTc concentration 
are found to be {cs, laTc) = 0 ng/ml}. The estimation 
of I/O response of genetic transistor based on experi- 
mental results is shown in Figure 2, with experimental 
details summarized in Additional file 1. Clearly, the I/O 
characteristics of genetic transistor can match the de- 
sired I/O response in a workable input range gi e [298, 
431] under the intrinsic fluctuations and environmental 
disturbances. 
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Figure 2 The amplifier design example of synthetic genetic 
transistor. For amplifier design example of synthetic genetic 
transistor, the prescribed amplification gain = -2 within the range 
G [298,431] is computed by (1 1) (red line). The most adequate 
promoter-RBS component C3 and aTc concentration laTc to fit the 
prescribed amplification gain are searched as {c^, laTc\ = {h,0 ng/ml} 
by minimizing JAiCsJaTc) in (16) from the corresponding promoter- 
RBS library Libret and concentration range of inducer lajc- The green 
points are the experimental results based on {C3, lard, and the error 
bars are the standard deviations. The green line is the estimation of 
I/O response of the synthetic genetic transistor based on 
experimental data, with the estimated amplification gain = -1.978. 
Obviously, the amplification gain of the designed genetic transistor 
could match the desired amplification gain quite well. 



Switch design example of synthetic genetic transistor 

Consider the switch design of the synthetic genetic tran- 
sistor. The switch design procedure is similar to the 
amplifier design procedure of a synthetic genetic transis- 
tor. Firstly, to obtain the complete I/O characteristics of 
switching, promoter-RBS components {ci, C2} = {/e, ^3} 
are selected to obtain the maximum input operation 
range. The dynamic model and the steady state model 
have been described in (2) and (3), respectively. The in- 
put operation range and desired I/O switch response are 
specified as follows: 
Input operation range: 



/lG[2*10-^10] (mM)^^iG[103,614] (a.u.) (14) 



and 



(3249.7-1,) 
l + (^i/150.4)^ 



(15) 



where denotes the low level of switching or basal level 
of promoter-RBS component C3. Note that the standard 
deviations of parameter fluctuations that are supposed 
to be tolerated in vivo and from environmental distur- 
bances are the same as in (12). Finally, GA is employed 
to search a set {c3,4rc} from corresponding libraries to 
minimize the following cost function: 



min Js{c3,hTc) 

C3GI^^^e^,/2G[2*10-^10] 

p614 

= min E 

cseLibTet M^[2*10-\10] J 103 



(yss (^3 JaTc.gi) -Jd {gl )fdg^ 
(16) 



Then, the most adequate promoter-RBS component 
from the corresponding library and aTc concentrations 
are found to be fej^rJ = {^3> 0 ng/ml}. The estimation 
of I/O response of genetic transistor based on experi- 
mental results is shown in Figure 3, with experimental 
details summarized in Additional file 1. Clearly, the 
switching I/O characteristics of synthetic genetic transis- 
tor can match the desired I/O response under the intrin- 
sic fluctuations and environmental disturbances. 

According to the above examples, the amplification or 
switching I/O characteristics of a synthetic genetic tran- 
sistor with different design specifications can be 
achieved by selecting the most adequate promoter-RBS 
component C3 and inducer concentration using the pro- 
posed library-based searching method. However, not just 
the promoter-RBS component C3 e Libxet can be selected 
to achieve the amplification or switching design specifica- 
tion of the synthetic genetic transistor circuit, but also 
other promoter-RBS components, Le,, Librae-, can be se- 
lected to achieve the desired I/O response. However, for 
various design specifications, more promoter-RBS libraries 
are needed to achieve these design specifications. 



2500 



2000 



1500 



^ 1000 
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100 150 200 250 300 350 400 450 500 550 600 650 
Input Fluorescence (a.u.) 

Figure 3 The switch design example of synthetic genetic 
transistor. For switch design example of synthetic genetic transistor, 
a desired I/O switch response is computed by (15) as shown in red 
line. The most adequate promoter-RBS component C3 and aTc 
concentration \ajc to fit the desired I/O switch response are searched 
as {C3, \aj^ = {Ti, 0 ng/ml} by minimizing J5(c3, \aj^ in (16) from the 
corresponding promoter-RBS library Ubjet and concentration range 
of inducer lajo The green points are the experimental results, and 
the error bars are the standard deviations. The green line is the 
estimated I/O switch response based on experimental data. 
Obviously, the I/O response of the designed genetic switch could 
match the desired I/O switch response quite well. 
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For the convenience of synthetic genetic transistor de- 
sign for synthetic biologists, one look-up table has been 
built for the various design specifications as shown in 
Table 1 via selecting adequate promoter-RBS compo- 
nents from the corresponding libraries and adequate in- 
ducer concentration to achieve the optimal matching in 
(9). Based on various amplification gains in some specific 
operation range, the synthetic genetic transistors can be 
designed by first checking the look-up table. In future, 
more promoter-RBS components and inducer concen- 
trations for different I/O characteristics of synthetic gen- 
etic transistors can be accumulated to build much larger 
look-up tables to match a lot of design specifications. 
From this look-up table, based on the desired design 
specifications, we can select the adequate promoter-RBS 
components and inducer concentrations to synthesize 
the genetic transistors with desired I/O responses. Thus, 
less time will be spent on the design procedure as a de- 
signer will be able to easily construct transistors with the 
desired I/O characteristics. 

Discussion 

One major aim of synthetic biology is to construct a 
gene circuit with the desired functionality of an organ- 
ism. Recently, promoter libraries and promoter-RBS li- 
braries have been built to simulate the in vivo behavior 
of a gene circuit [30,38,39]. By identifying the kinetic 
strengths of promoter-RBS components, the protein 
expressions in the gene circuit can be estimated and 
predicted. However, in the process of constructing 
promoter-RBS library, the identified kinetic parameters 
in the promoter-RBS library can be affected by several 
conditions, including the medium, copy number of plas- 
mid, terminator and so on. Therefore, for the extensive 
application of promoter-RBS libraries, the construction 
conditions of promoter-RBS libraries need to be unified 



Table 1 The look-up table with different gain 
specifications for synthetic genetic transistors 

Amplifier gain specifications of synthetic genetic transistor 





Input range 
(a.u.) 


Gain 


Promoter-RBS 
component 


Inducer 
concentration 


Libjet 


120-180 


-10.00 


T3 


0 ng/ml 




150-225 


-7.50 


T3 


1 ng/ml 




260 ~ 460 


-2.00 


h 


0 ng/ml 




400-550 


-1.00 


T3 


1 ng/ml 




460 ~ 560 


-0.75 


T3 


0 ng/ml 




575 -620 


-0.50 


T3 


1 ng/ml 


Libiac 


40-140 


-0.15 


^1 


0 mM 




40-140 


-2.50 


^3 


0 mM 



Given the desired amplification gains and their input signal ranges, we could 
select adequate promoter-RBS component and inducer concentration from the 
table to achieve the minimum matching error in (9). 



and standardized. This will allow standardized promoter- 
RBS libraries, similar to electronic component libraries, 
which can be easily used and expanded by other gene 
circuit designers. 

In this study, by the promoter-RBS libraries we estab- 
lished, a genetic transistor has been constructed and im- 
plemented. Additionally, the synthetic genetic transistor 
can perform amplification and switching like an elec- 
tronic transistor according to its I/O characteristics. 
The I/O characteristics of the synthetic genetic transis- 
tor circuit are simulated by a mathematic model with 
random parameter fluctuation to guarantee the robust- 
ness of the design in vivo. The design specification of 
amplification or switching in the genetic transistor can 
be achieved by the library- searching method using GA. 
By optimally matching the desired I/O response 
of amplification or switching, the most adequate set of 
promoter-RBS component and inducer concentration 
{c3,/2} can be selected to construct a genetic transistor 
with the desired design specifications. The library- 
searching method using GA is introduced to reduce 
the number of trial-and-error attempts, as well as the 
searching time in libraries when the libraries have a 
large number of components. Furthermore, for the 
convenience of synthetic genetic transistor design for 
synthetic biologists, one look-up table has been built 
for the various design specifications as shown in 
Table 1. From this look-up table, based on the desired 
design specifications, we can select the adequate 
promoter-RBS components and inducer concentrations 
to synthesize the genetic transistors with desired I/O 
responses. Thus, less time will be spent on the design 
procedure as a designer will be able to easily construct 
transistors with the desired I/O characteristics. 

For applications of the genetic transistor, the various 
biological components need to be characterized and 
standardized. By using characterized and standardized 
genetic components, the design specification of a genetic 
transistor can be set and the look-up tables can be 
used to support the genetic circuit design. The genetic 
transistor described here has a number of potential 
applications. The amplifier can be used to amplify the 
oscillation signal reversely and linearly. Based on the 
designed oscillatory genetic circuits [7,11-13,40-42] in 
oscillatory metabolic pathways [43-45], an adequate gen- 
etic transistor selected from the look-up tables according 
to the oscillation range and desired amplification gain 
can be inserted into these circuits directly to amplify the 
oscillatory signal. In this way, the original genetic cir- 
cuits do not need to be redesigned. On the other hand, 
the switch can be used to detect some signals and act 
like a detector or biosensor [25,29,46,47]. When the in- 
put signal changes, the output signal will switch to the 
other state and make the downstream circuit respond to 
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the signal change. Also, the switch of a genetic transistor 
can work as logic gates as in an electronic transistor 
[16-18]. With different combinations of genetic transis- 
tors, different logic gates can be constructed. 

Conclusions 

In this study, three kinds of libraries, i.e., a constitutive 
promoter-RBS library, repressor-regulated promoter- 
RBS library and activator-regulated promoter-RBS li- 
brary, were established for constructing synthetic gene 
circuits with the desired transistor amplification or 
switching function. The amplification gain and switching 
level of a genetic transistor could be calibrated by select- 
ing adequate promoter-RBS components and inducer 
concentrations from the corresponding libraries. For the 
measurement of I/O response, we could embed an add- 
itional repressor-regulated promoter-RBS component 
with reporter protein at the input terminal to measure 
the input signal while replacing the output protein with 
a reporter protein to measure the output signal. Further, 
for the convenience of input regulation, an external cir- 
cuit was constructed to control the input signal using 
the concentration of inducer. Based on the desired I/O 
response in relation to amplification or switching of a 
genetic transistor, the GA-based searching algorithm 
was introduced to search for the most appropriate set of 
promoter-RBS components and inducer concentration 
from the corresponding promoter-RBS libraries to achieve 
prescribed I/O characteristics in the genetic transistor. In 
the simulation results for this study, we demonstrated that 
the genetic transistor designed here has the prescribed 
function of amplification or switching. By the library- 
searching method using GA, different design specifications 
of amplifier or switch could be achieved the most appropri- 
ate set of promoter-RBS components from the correspond- 
ing promoter-RBS libraries and inducer concentration 
within a feasible region. Finally, a look-up table was built 
for genetic transistor design with different genetic transistor 
design specifications. Using this table, we could easily select 
an adequate set of promoter-RBS components and inducer 
concentration to construct the desired genetic transistor. 
This innovation saves much time in trial and error attempts 
in the iterative experimental procedure. 

Additional file 



Additional file 1: Supplementary Appendix. 
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